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Abstract 

We study the vahdity of the large- A?^ equivalence between four-dimensional SU (N) lattice gauge 
theory and its momentum quenched version — the Quenched Eguchi-Kawai (QEK) model. We find 
that the assumptions needed for the proofs of equivalence do not automatically follow from the 
quenching prescription. We use weak-coupling arguments to show that large-A^ equivalence is 
in fact likely to break down in the QEK model, and that this is due to dynamically generated 
correlations between different Euclidean components of the gauge fields. We then use Monte- 
Carlo simulations at intermediate couplings with 20 < N < 200 to provide strong evidence for 
the presence of these correlations and for the consequent breakdown of reduction. This evidence 
includes a large discrepancy between the transition coupling of the "bulk" transition in lattice gauge 
theories and the coupling at which the QEK model goes through a strongly first-order transition. 
To accurately measure this discrepancy we adapt the recently introduced Wang-Landau algorithm 
to gauge theories. 



I. INTRODUCTION 



QCD simplifies in the 't Hooft limit of a large number of colors, yet differs from the 
physical case, = 3, by only 0{1/N) corrections [D [2]- As a result, it has been a long- 
standing goal to understand the properties of QCD at large N . This goal has become even 
more interesting as string theory developments, based on gauge/gravity duality, attempt 
to move towards predictions for QCD-like theories at infinite [3J. Considerable progress 
towards this goal has been made using numerical lattice computations (as reviewed, for 
example, in Refs. [H [3]). For example, using conventional large volume simulations, and 
extrapolating from A^ = 3 — 16, precision results for several quantities in the large- A^ limit 
have been obtained. In addition, many quantities have been found to depend only weakly 
on A^. 

In this paper we reconsider an alternative to conventional large volume simulations, 
namely the use of large-A^ volume reduction, in which the infinite space-time volume is 
collapsed to a single point, with space-time degrees of freedom repackaged into the 0(A^^) 
color degrees of freedom. This allows one, in principle, to trade two large parameters, the 
volume and A^^, for a single large parameter, and consequently to consider much large values 
of N. 

The idea of large-A^ reduction for lattice gauge theories was first proposed by Eguchi 
and Kawai in Ref. [6j. They defined the "reduced" theory to be a lattice gauge theory 
on a single site and showed that, under certain assumptions, Wilson loops in the reduced 
theory and in the d-dimensional infinite-lattice gauge theory acquire the same expectation 
values in the large- A^ limit. The two main assumptions in [6] are that expectation values 



of products of single-trace operators factorize at large- A^ (see Eq. (2.9) below) and that 



the vacuum is symmetric under the {Z^Y center transformations applied to the model's 



SU{N) "link" matrices [/^ (see Eq. (2.5))^. It was quickly realized, based on weak coupling 
arguments and numerical results, that the second assumption does not hold for > 2 in the 
continuum [71 El [9], and various ideas for solving this problem were suggested. The first, 
which we focus on here, was of quenching the eigenvalues of the link matrices by forcing them 



^ We consider the SU{N), rather than the U{N), theory in this paper. The two theories become equivalent 
as iV ^ oo, but the extra phase in the U {N) theory, which decouples from the dynamics at any N, obscures 
the underlying mechanisms we discuss in later sections. 
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to have a Zjy invariant distribution. This means that, by construction, all order parameters 
for the breakdown of the {Zj^iY symmetry, such as (trf/^), vanish. 

This "quenched Eguchi-Kawai" (QEK) model was proposed in Ref. |7j (see also [101), and 
first numerical results indicated that it did indeed solve the problem of unwanted symmetry 
breaking [H], [121 lEl ttH [iSl [16]. Intensive analytic investigation of the QEK model ensued. 
(For example, see the papers [171 [TU [191 EOl [211 [22], which are relevant for us here.) For 
further discussion and references we refer to the reviews in Refs. [231 1211 [25] . 

Studies of the QEK model tailed off, partly because of the emergence of the apparently 
more promising "twisted Eguchi Kawai" (TEK) model [2S]- We do not discuss this model 
here, but note only that despite the early literature, recent extensive simulations find evi- 
dence that large- iV reduction fails in the TEK model [271 EHl [29] . 

Another line of development, initiated in Ref. [30] and extended in a series of papers 
summarized in [5], has been much more successful. This is the idea of partial reduction. 
Here one reduces not to a single-site but to a lattice of size L'^, with L > L^^^^ ^ 1 fm. 
As long as the dimensions exceed this minimum, the {Z^Y center symmetry is seen to be 
unbroken, so reduction holds and one obtains volume-independent results if is sufficiently 
large. Detailed studies have demonstrated partial reduction, and provided results for gauge 
and fermionic quantities for A^ < 50 [5j . 

In this paper we return to the original QEK model, and study whether reduction holds 
there as well, as was indicated by the early works [TU |T21 HSl HH UHl US] • We are motivated 
to do so not just as a tool to study very high values of A^, but also for the following 
reason. It has recently been argued that the {Z^Y breaking does not occur in any volume 
if the fermions are in the adjoint representation of the gauge group, and that in such a 
theory, complete large- A^ reduction to a single site should hold [31]. This theory is of 
phenomenological interest because, through the orientifold large- A^ equivalence of adjoint 
and two-index antisymmetric tensor fermions, it differs from physical QCD by corrections of 
0(1/A^) [32j. Studying this theory numerically at any volume, however, is expensive since 
it requires simulating dynamical fermions.^ Thus, before plunging in to such a study, we 
chose to get experience with single-site models, and the QEK model is a natural candidate. 

^ In contrast to the 't Hooft limit where the fermions are in the fundamental representation and affect the 
dynamics at 0(1 /N), in this theory they enter at 0(1) and thus dynamical simulations are necessary. 
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The outline of paper is as follows. In Section |TT] we review the features of the original 
Eguchi-Kawai (EK) model which are relevant for our study, such as the failure of reduction 



in this model, and its relation to (Z^) breaking. In Section III we discuss the QEK model 
in some detail. We begin by introducing the quenching prescription, and by presenting 
the symmetries of the model. We then review the theoretical arguments for the validity 



of the large-iV equivalence between the QEK and large-A^ QCD. In Section IV we present 
weak-coupling analytic considerations that question these arguments and we conclude that, 
similarly to the EK model, it is likely that reduction fails in the continuum limit. In Section fV] 
we present an extensive numerical lattice study of the QEK and find evidence to support 



our claims of Section |IV| We conclude in Section |VI] with remarks on the implications of our 
findings for other large- reduced models, such as those in |33] . 



II. REVIEW OF THE ORIGINAL EGUCHI-KAWAI MODEL 

We begin with a brief review of the original EK model, so as to set notation and provide 
background for our observations. The EK model is a d-dimensional lattice gauge theory 
restricted to a single site whose partition function is 

Zek = J DUexpiSEK), (2.1) 
Sek = NbJ2 2Re Tr {U^UMIuD . (2.2) 

Here b = 1/X, with A = g^N the 't Hooft coupling, and the integration measure is the Haar 
measure on SU{N). Aside from the "reduced" gauge symmetry, 

V/x : u^^nu^n^ with nesu{N), (2.3) 

the model is also symmetric under center transformations applied independently to the d 
link matrices 

Uf, -> Uf, z""^ with z = e^^'^ and G Z^v , (2.4) 
and under the d reflections 

-P,: U.^Ul (2.5) 

A Wilson loop on the original lattice that is defined by the path C : (x, a; -|- /i, . . . , a; — 
z> — p, a; — z>, x), and that is given by 

Wc = ^ti f/a-,A^a;+A,!> ■ ■ ■ U^~u-p,pUx-0,u , (2.6) 
4 



is mapped in the reduced model to 

^c" = ^ tr U^U, ■ ■ ■ UpU, . (2.7) 

The essence of reduction is that the expectation values of Wc in the gauge theory and of 
j^reduced ^^^q rcduccd modcl are the same at large- iV:^ 

(M^c)gauge theory = (1^5'""^'=') reduced + 0{1/N^). (2.8) 



To obtain Eq. (2.8) one derives the Dyson-Schwinger equations of Wc and Wq'^^'^'^'^ and finds 
that they are the same, up to additional source terms. These source terms vanish provided 
that the following two conditions are satisfied 

{Wc,Wc,) reduced = (Wc,) reduced (Wc,) redneed + O {1 / N^) , (2.9) 
(W^open) reduced = 0. (2-10) 



Here Eq. (2.9) must hold for all contours Ci and C2, while in Eq. (2.10) W^open denotes any 



reduced Wilson loop whose path is a mapping of an open path in the gauge theory, e.g. 

Copen-ix,x + fi,...,y-i',y) ; y^x. (2.11) 



This means that, in at least one direction, the corresponding H^onen'^^'^ of Eq. (2.7) has the 



open 

the number of U's less the number of W's different from zero."^ 



The factorization required in Eq. (2.9) is expected to hold in large- theories, and realizes 
the idea of the "Master field" |2]. Condition (2.10) holds if the vacuum is {Z^Y symmetric. 



since under this symmetry lyopen acquires a phase. In fact, as already noted in Section |T], 
this symmetry is spontaneously broken in the four- dimensional EK model, at weak-coupling. 
This is shown by perturbative calculations in the weak-coupling (6 00) limit [71 [H |3l] , 
and has been demonstrated by lattice simulations to hold for b > 0.19 — 0.30 [TJ |9]. Thus 



Eq. (2.10) does not hold, invalidating reduction. 

It is useful for our subsequent discussion to briefiy recall the perturbative calculation of 
Refs. [71 [8] . One first writes the link matrices in polar form 

U, = ViA,V,, (2.12) 



^ See, however, the discussion of connected correlation functions in [31] . 

^ More precisely, this difference need only be zero modulo N, a subtlety that will not play any role in our 
considerations. 
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where G SU{N) and is a diagonal matrix containing the eigenvalues, 



= diag[exp(ipj,), . . . ,exp{ip^)] ; p"^ G [0, 27r) . 



(2.13) 



For Ufj, G SU{N), the are constrained to satisfy^ 



N 



E P)l = mod 



(2.14) 



a=l 



in each direction — a constraint that we keep implicit in the following formulae for the sake 



of clarity. Using Eq. (2.12) one can show that the partition function Eq. (2.1) becomes 



Zek = JY[j^ ^\p) J DV exp SEK^j\{{^ exp [-Fek{v)] 



(2.15) 



where 



M a<b 



(2.16) 



is the Vandermonde determinant, and DV is the Haar measure on SU{N). For all values 
of p^, the action Sek is minimized when = 1 for all fi (up to gauge transformations). 
Expanding around unity, and assuming non-degenerate p, i.e. 7^ if a 7^ b, one finds 
at leading order in the weak-coupling expansion that the free energy is [71 [S] 

^6 



FEK(p)'^(ci-2)Elog 



a<b 



E 



sm 



Pl-P], 



(2.17) 



For d > 2, Fek{p) is minimized when, for each /x, the for all a become equal. ^ This 
implies that, for SU{N), the theory has A^'^ "vacua", with 11^ ~ e*^!^! ~ which are 



transformed into each other by the center transformations (2.4). 



Such a clustering of the eigenvalues appears to indicate spontaneous breakdown of the cen- 



ter symmetry and so to invalidate Eq. (2.10). To establish spontaneous symmetry breaking 



(SSB), however, one needs to know whether fluctuations about each vacuum are sufficient 
to restore the symmetry (as happens for < 2 in infinite volume statistical mechanical 
systems). In other words, does the free-energy barrier between the different vacua become 



^ We use a as a color index, and denote the lattice spacing by aiat- 

^ As noted in Ref. 7 , the weak-coupling calculation of F^k breaks down for the degenerate eigenvalues 
that are picked out by minimizing -Fek- This is due to the presence of zero modes. One can extend the 
calculation to include the effects of these modes and the conclusions are unchanged in the weak-coupling 
limit [5H. At moderate values of b, however, only numerical calculations are reliable, and these [SHU] are 
consistent with the weak-coupling picture of ^. 



6 



infinitely high as — > oo (implying symmetry breaking) or not (implying symmetry restora- 



tion)? Also, the calculation leading to Eq. (2.17) is valid only when b oo, and it is possible 
that the symmetry is restored for moderate values of b. In fact, as noted above, numerical 
simulations imply that the symmetry is indeed broken once b becomes moderately large. 
We conclude this section with a general remark concerning the possible ways that re- 



duction can fail. Of the two key conditions, Eqs. (2.9)-(2.10), it is often considered to be 
the second that is crucial. Thus the validity of Eguchi-Kawai reduction has become almost 
synonymous with the absence of spontaneous breaking of the center symmetry. In this 
scenario there are multiple vacua, connected by symmetry transformations, around each of 
which the fluctuations are ~ 1/N'^. We wish to emphasize, however, that this is not the 



only possibility. What is required for reduction to hold is the combination of Eq. (2.9) and 



Eq. (2.10), and it is also possible that the first of these can fail. This happens if there are 
multiple would-be symmetry-breaking ground states yet fiuctuations lead to motion between 
all these states even when ^ oo. (This possibility has been already mentioned in the 
previous paragraph.) In an infinite volume theory this corresponds to a breakdown of clus- 
ter decomposition. These two scenarios should be compared to that with a single vacuum 
obeying cluster decomposition, in which case both relations hold and reduction is valid. 

The distinction between the failure of reduction due to a breakdown of cluster decomposi- 
tion and due to SSB is important below, so we illustrate it with the following simple example. 
For an SU{N) EK theory, as noted above, there is SSB for sufficiently weak coupling, with 
A^'^ vacua. If we change the gauge group to U{N), however, the free energy landscape has 
a "Mexican-hat" form, and the A^'^ vacua become part of a continuous degenerate manifold 
connected by the four U{1) phases (one per direction). The path integral over these phases 



causes expectation values of open loops to vanish and Eq. (2.10) holds. This does not, how- 



ever, mean that reduction holds, because the first of the required conditions, eq. (2.9), is 



not satisfied. The failure of (2.9) can be seen by comparing its two sides for the case where 
Ci and C2 are both open loops that, when combined, form a closed loops (for example see 
Fig. [T|. For this choice the U{1) transformations multiply Wc^ and by opposite phases 



and so the l.h.s. of Eq. (2.9) is independent of these U{1) phases and has an 0(1) value. In 
contrast, the r.h.s. is of 0(1/A^^) because the open-loop expectation values do vanish. The 
failure of factorization is perhaps surprising, but occurs because the vacuum is not unique 
and because the U{1) degrees of freedom have unconstrained fiuctuations. This is also an 
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example of why it is simpler to analyze the SU{N) theory. 



III. REVIEW OF THE QUENCHED EGUCHI-KAWAI MODEL 

Quenching attempts to avoid the spontaneous breaking of the {Z^Y symmetry by forcing 
the eigenvalues not to cluster. In this section we recall the quenching prescription and 
describe its relation to infinite-volume large-A^ QCD. 



A. The quenching prescription 

The prescription consists of first calculating expectation values for a fixed set of the 
eigenvalues (labeled collectively by "p" and distinguished from the usual expectation 
values by a subscript), 

{0{U))^ = Z{p)-' [ l[DV^ e^Q--(p) 0{U) . (3.1) 

Here S'qek(p) is simply ^ek with the having the form ( |2.12[ ), i.e. 

^qek(p) = iV& E 2Re Tr {(ViA.V,) (V^AM iViA%) (KIA+K) ) , (3.2) 

Z[p) is the partition function for fixed p, 

Zip) ^ [ l[DV,expiSciEK{p)) , (3.3) 

and the take values in SU{N). The second part of the prescription is to average expec- 
tation values over the choices of p: 

(0(W))qek = [dp (0(f/)), = r n ^ Pip) {Om, . (3.4) 

•' "'0 ij.,a 

Here p{p) is a positive weight function (with integral normalized to unity) which is invariant 
under the {Z^Y shifts Pu ^ Pu^ 27rA;^/A^, and dense in the space of the as — > oo. For 



the SU{N) theory, it should also incorporate the constraints (2.14). We discuss particular 
choices of p{p) below. 

To understand the significance of quenching, it is useful to write expectation values in 
the original EK model in terms of the p-dependent {0{U))p: 

8 



Comparing Eq. (3.5) and Eq. (3.4) we see that quenching changes the measure of the integral 
over the p in such a way as to replace the non-uniform weighting A'^{p)Z{p) = exp[— Fek(p)], 
which was the cause of the clustering of eigenvalues |TU], with the uniform weighting p{p). 



B. Symmetries in the quenched model 

Since quenching has divided the original dynamical degrees of freedom (the U^) into 
the dynamical and the quenched it is important to understand how the symmetries 



Eqs. (2.3)-(2.5) are realized. The gauge transformations Eq. (2.3) can be chosen to act only 



on the V^: 



(3.6) 



By contrast, the center and reflection symmetries must, in general, be realized by transform- 



ing the eigenvalues: center transformations (2.4) become 

27m 



p^^mod(p^ + ^,27r 



while the reflections (2.5) become 



(3.7) 



(3.8) 



There are, however, special choices of p{p) for which one can realize Eqs.( 2.4 2.5) by trans- 



formations on the Vn alone, and we discuss these in Section III C below. 



Quenching solves the problem of the dynamical clustering of eigenvalues, and leads to 
the desired vanishing of the expectation values of open Wilson loops such as tr (Ufj,). To see 



this, note that the center transformation Eq. (3.7) performs a "clock rotation" of the and 



thus multiplies (tr(f/^))p by z"''. This is as in the unquenched EK model, but now the 
are forced to have a Zf^ symmetric distribution. Thus the integration in Eq. (3.4) leads to 



QEK 



dp {triU,))^ 



0. 



(3.9) 



In fact, the quenched expectation value of any open loop will vanish due to the (now enforced) 
average over the center transformations. 

C. Large-A^ reduction in the QEK model 



As in the EK model, one can derive the Dyson-Schwinger equations for Wilson loops in 
the QEK model. They too include unwanted source terms. In the QEK model, some of 



A 



FIG. 1: Illustration of the extra source terms, eq. (3.10), in the Dyson-Schwinger equations for 



a closed Wilson loop in the QEK model. The dashed line represents Wopen while the solid line 
represents W^pen- 



these have the form 



/ dp (PVopen 



open/ 



(3.10) 



where PVopen and W'^ are open Wilson loops (as defined in Section II) that, when joined, 



form a closed loop, as illustrated in Fig. [T| Such terms must all vanish for reduction to hold. 
The argument that they do vanish proceeds in two steps [171 [TBI [19] : 



/ dp (l^open 
/ dp (I^open), (W^^ 



/ dp (I^open), (l^^p,,)^+ 0(l/iV2) (3.11) 
dp (W^open), / dp' (W^^pen) , + 0(l/iV) . (3.12) 



The first step is large- factorization for a fixed set of p, valid to all orders in perturbation 
theory. The second step, which one might call "quenched factorization", is special to the 
quenched theory. If it holds, then, due to the vanishing of quenched expectation values of 



open loops [as in Eq. (3.9)] the extra terms (3.10) in the loop equations do vanish in the 
large- limit. 



We will argue in subsequent sections that the combination of eqs. (3.11) and (3.12) does 



not hold, most likely due to a failure of the latter equation. In order to understand what 
fails, we must first describe the argument for the correctness of these steps, and indeed why 
the quenching prescription is expected to reproduce the large- dynamics of SU (N) lattice 
gauge theories. 



^ Here we note that the derivation of the Dyson-Schwinger equations in the QEK model is different than in 



the EK model, and in addition to the terms of the form of Eq. (3.10 1, there are other source terms which 
have a similar, but more complicated form |18j . Nevertheless, the analysis we perform in this section 
holds for these terms as well. 
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The approach of Refs. |71 [TTl [HI [191 [2Q] provides an intuitive explanation of why large- 
quenched reduction works for a wide class of theories, albeit within perturbation theory. The 
idea relies on the fact that, in large-A^ perturbation theory, only planar diagrams survive. 
In 't Hooft's double-line notation each gluon line is replaced by two oppositely pointing 
lines that carry two indices (a, b) with a,b E [1, A^]. In the reduced theory there is, initially, 
no momentum associated with a gluon "propagator". The key point is then to associate a 
(i-dimensional lattice momentum with each of the indices, 

a ^ p",; p;:G[0,27r), (3.13) 

and to assign to the gluon with color indices (a, b) the difference in the momenta associated 
with the two indices: 

g^" ^ mod (p;: - , 27r) . (3.14) 

It is easy to check that since one lets all take all values in [0, 27r), the momenta q^^^ of all 
gluons in any planar diagram will independently take values in the Brillouin zone, and will 
obey momentum conservation (modulo 27r) at the vertices. This is impossible for non-planar 
diagrams, where some of the gluons carry the indices (a, a) and so have q°f' = 0. 



The identification of Eq. (3.13) is necessary in order to embed space-time (or rather the 
first Brilluoin Zone of its momentum space) in color space, but it is only the first step. The 
next is to choose the action of the single-site model in such a way that the actual value 
of planar diagrams will be the same as that in the full gauge theory. In particular, they 
choose it such that the (a, b) matrix element of the gluon propagator takes (in an appropriate 
gauge) the usual "1/g^" form (or, more precisely, its lattice version) with q indeed being 
the difference ]f — p^. Vertices are similarly reproduced. In this way, for a given choice of 
color indices, one obtains the correct integrand of the corresponding infinite- volume large- iV 
Feynman diagram. 

The application of the quenching prescription involves an important subtlety, which re- 
quires that additional constraints be placed on the single-site fields [HI [20]. The end result is 



that one arrives at precisely the QEK model, with the quenching prescription of Eqs. (3.1 )- 



(3.4). The quantities p° are now viewed as (dimensionless) loop momenta. For example, 



perturbing around the classical vacuum of Ki = 1, yields the following two-point function 
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for the fluctuating matrices 



(3.15) 



E.sin2(2^) 

which is the standard lattice result for the gluon propagator.^ 

The final step in the quenching prescription of Ref. [3 [TSl [121 ED] is to average over the 



momenta, as in eq. (3.4). The weight function should be manifestly {Zj^jY invariant, and, as 
N ^ oo, force the momentum-components in each direction to densely cover the Brilluoin 
zone. An example for such a measure is that suggested in : 



PVdM(p) = A2(p) 



(3.16) 



where the Vandermonde determinant is defined in Eq. (2.16). At large- the function 



PYduip) forces the momenta to lie as far apart from each other as possible, and combined 



with the SU{N) constraints (2.14), this requires the to be a permutation of the "clock" 

A^ + 1- 



values, 



pa 



27r 



a — 



a e [1,N] 



(3.17) 



A^ V"^ 2 

The integral over p then amounts to an average over permutations, independently for each 



direction. Since the momenta (3.17) are uniformly distributed, this gives a discrete ap- 



proximation to the integration Eq. (3.4) over the infinite- volume momentum space. Thus 



as A^ — s> oo one obtains, order by order in perturbation theory, the correct infinite-volume 
result for each Feynman diagram. 

We can now give the argument of Refs. [TSl HH] for the crucial relation (3.12). Imagine 



evaluating the l.h.s. of Eq. (3.12) in perturbation theory, and focus on the contribution from 



planar diagrams with (L+M)-gluon loops, with L loops coming from the expansion of PVopen, 
and the remaining M loops from the expansion of W^open- We can write this contribution as 



dp fiPai,Pa2,---,PaJgiPb^,Pb2,---,Pb 

" m .nn n. t 



(3.18) 



bl,b2,.-.,bM 



Here / and g denote the integrands of the planar diagrams contributing to the two Wilson 
loops, and pa^ and pb^ are the momenta that flow in these diagrams. If {L + M) <^ A^ 



^ As explained in [HIEO], the expansion used to obtain Eq. (3.15) is V^A^V^Aj^ = e*""^''. 
^ Note that here the Brilluoin zone is, at large-iV, [— tt, tt)'^, instead of [0, 2ttY, but this change is irrelevant 
in our discussion. 
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then a generic term in the double sum has all indices different, and for each such term the 
(normalized) integral over p factorizes in the large- limit into the two integrals 

j dp f{Par,Pa2, ■ ■ ■ ,PaL) giPbi,Pb2, ■ ■ ■ ,PbM) = J dp f (Pa^, Pa2, ■ ■ ■ , PaJ J dq g{qb„ Qb^, ■ ■ ■ , qbu'. 

(3.19) 

Here, for clarity, we changed the dummy integration variables pb^ to qb^. The last step is 



to sum Eq. (3.19) over all possible values of the indices and bi, which gives the L-loop 



contribution to (PVopen) multiplied by the M-loop contribution to (W^opcn)- If this step is 



correct, we obtain the r.h.s. of Eq. (3.12). 



This last step is only approximately correct for the following reason: Eq. (3.19) holds 
only if the indices ai,,..,L are all different from the indices while in performing the 

final sum we ignore this restriction. At large- A^, however, the effect of this "negligence" is 
only of 0(1/A^) — the fraction of terms with equal and bi indices. As a result one finds 
that to all orders in planar perturbation theory, quenched reduction holds at large- A^. The 
fact that one ignores 0(1/A^) terms here also explains why the quenched model has 0(1/A^) 
rather than 0(1/A^^) corrections. 

D. Alternative choices of p{p) 

The formulation of the QEK model provided by the approach of Refs. [TSl HO] shows that 
there is considerable freedom in choosing the weight function p{p). The choice must simply 
turn the integrands of Feynman diagrams into their integrals as A^ oo. The simplest 
choice is to use a uniform distribution: Puniform(p) = 1 for all A^. In practice, this can be 
implemented by Monte-Carlo — drawing p randomly from a uniform distribution. Another 
choice, which we call Pciock(p); is to take the momenta to be a permutation of the clock values 



of eq. (3.17), even for finite A^, and then average over permutations. This corresponds to 
setting 

pcioM = (tie) {um - p^"^^"^)} ' (3-20) 

where the (J^{a) are permutations of the color indices. With this choice, which we use ex- 
tensively below, the center-symmetry-breaking order parameters tr UJ^ (for p ^ mod A^) 
vanish prior to the average over permutations. This choice also has a simple physical inter- 
pretation. The discrete are those that one would obtain if one had a lattice with A^ sites 
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in each direction, with periodic/antiperiodic boundary conditions for odd/even. After 
averaging over permutations one obtains the result of the Feynman diagram on a lattice of 
physical volume (A^aiat)'', where aiat is the lattice spacing. 

It is further argued in Ref. [HI [19] that the integral over momenta (or sum over permu- 



tations) becomes unnecessary as — > oo. This is because Eq. (3.14) implies that the sum 
over the color indices in each color loop becomes an integral over the corresponding Brillouin 
zone. Thus a single choice of randomly chosen momenta, or a single set of randomly chosen 
permutations of clock momenta, should, in principle, be sufficient. In practice, for finite N, 
it may be preferable to include an average over such choices. 

The final choice of p{p) we consider is that suggested in [18] and analyzed by Bars in [35] . 
It applies only when = K'^, where K is an integer. In this case p{p) is a product of delta- 
functions such that each value of the color index is associated with a different d-dimensional 
momentum lying on a K'^ latticization of the Brillouin zone (BZ). In four dimensions one 
has 



{Pl, P2, Pi Pi} = 

27r r 

— < mod (a — 1, K) , mod 
K I 



K 



K I , mod 



K I , mod 



K ,(3.21) 



where [x\ indicates the integer part. We shall henceforth denote this choice by "BZ". As 
an example, consider the two-dimensional case with = 16 (so K = 4). We divide the 
Brillouin zone into sixteen boxes, and set the sixteen momenta p"" to lie at their centers, as 
shown in Fig. [2j 

(0,2ji) (2;i,2jt) 
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FIG. 2: A two-dimensional example of the embedding suggested in Refs. |18|I35| of the color indices 
in the Brilluoin zone for = 16. 

For p{p) = Pbz{p), the momenta are those for a physical volume {Kaig^tY = ^('^fat); 
which is much smaller than that obtained using the clock momenta. The advantage of using 
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P = Pbz is that one obtains a uniform distribution over the Brillouin zone from a single set 
of momenta, aheady at finite A^. 

The general discussion of the realization of center and reflection symmetries given in 



Section IIIB must be modified for pciock and pez- These symmetries can now be realized 



by a transformation on the V^. Consider first the clock momenta. There is then an SU{N) 
matrix S, such that, when V^j — > S'V^, 

= V^A^V,, ^ V^S^A^SV^ = Vi{A^z^)V^ = U^z^ . (3.22) 

This is possible because, first, up to an overall phase, the elements of are a permutation 
of the N'th roots of unity (so that multiplication by corresponds to a permutation of these 
elements) and second the eigenvalues of can be arbitrarily permuted by conjugation by 
SU{N) matrices (as will be discussed in detail in the subsequent section). The reflection 
transformations also correspond to permutations of the eigenvalues, and can be accomplished 
by different choices of S. 

The situation is slightly different for pez- Here the quenched theory only realizes a {ZkY 
subgroup of the center symmetry, because only such transformations correspond to a per- 
mutation of the momenta. The reflection transformations are also realized by permutations. 



IV. BREAKDOWN OF QUENCHED REDUCTION - ANALYTIC CONSIDERA- 
TIONS 

As explained in the previous section, the validity of quenching is predicated on the mo- 
menta being fixed by hand, independently in each direction, and then integrated over with 
a suitable weight function. This is possible in perturbation theory. When one does a 
non-perturbative calculation, however, the values of p are not completely fixed, and their 
distribution is thus determined in part by dynamics. In other words, they are incompletely 
quenched. We argue in this section that, at least in the weak coupling limit b ^ oo, the 
dynamics is likely to choose a ground state in which this incomplete quenching invalidates 



reduction in perturbation theory, and also invalidates the key relations (3.11 3.12). If this 
persists beyond perturbation theory, then reduction fails. Our numerical results, obtained 
for finite b, suggest that this is indeed what occurs. 
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The key observation is simply stated. The "fixed" momenta can be dynamically per- 
muted, independently in each direction, by fiuctuations in the V^. These innocent-sounding 
permutations lead, in general, to a different free energy, and the dynamics chooses the per- 
mutation(s) with the lowest free energy. The p's that one puts in by hand are different, in 
general, from those chosen by the dynamics, and the latter are not uniformly distributed in 
the Brillouin Zone. Thus, the sum over color indices does not lead to a uniform integration 
over the Brillouin zone, and the agreement with infinite- volume perturbation theory fails. 

The presence of permutations in the dynamics has long been known, and was stressed 
particularly in Refs. [22] and (SO]- To our knowledge, however, the implications for the 
validity of reduction have not previously been noted. 



A. "Momentum locking" at weak coupling 

Permutations are generated by transpositions, in which for one pair (ab) and 

one choice of fi. Transpositions can be accomplished, for example, by multiplying from 
the left by the SU{N) matrices V^"'''\(f)), which are the identity aside from 

f\/(«^)U)) =6''^"'. (4.1) 

V ^^'J {ab) block ^ ' 

As runs from to 7r/2, V^K^V^ changes to V^^A'^V^, where A^ differs from A^ in having 
the a'th and fe'th momenta permuted. Since it is possible to reach any permutation by a 
sequence of transpositions, an ergodic simulation will pass through all possible permutations 
of the input momenta. The question then is which of these permutations has the smallest 
free energy. Only if they are equally likely will the quenched model work as desired. 

We can calculate the relative free energies in the weak-coupling limit. First we describe 
the energy (i.e. action) "landscape". The minimum energy states, after appropriate gauge 
fixing, have = 1 and the momenta in any permutation of their input values. This is 
because the plaquette is unity for any choice of diagonal A's. If the momenta are non- 
degenerate, there are in fact (A^!)'^~^ different "vacua" (one factor of A^! being removed 



using a gauge transformation Eq. (2.3) to keep A^=i in its input order). These vacua are 



connected by the ^^"^•'(0), with the energy barrier (at = 7r/4) being [2T| [22] 

- ASqek = SNb sm\Apf/2) ^ sm\Apt'/2) . (4.2) 
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Here the transposition is being done on A^, and Ap^^ = P^i ~ p'ji- Generically, all the Ap 
are of 0(1), and the barrier height then grows as A^. There will always be some {ab) pairs, 
however, that have Ap = 0{1/N), and for these the energy barrier vanishes with increasing 
N. Nevertheless, for fixed N, as b oo, these barriers too become infinitely high. 

Thus, in the weak-coupling limit, and assuming non-degenerate momenta, one can treat 
the system as a collection of independent vacua with fluctuating around unity in each. 
At leading order the free energy in each vacuum is, up to an irrelevant constant, [71 18] 

-\nZ{p) = FeM + F^ip) , (4.3) 
Fek{p) = {d-2)Y: log (e sm'iApfm] (4.4) 

a<b V M / 

F2(p) = -\nA\p). (4.5) 



This is just a repetition of the result already quoted in Eq. (2.17), taking into account the 
difference in the definitions of Z{p) and Zek{p)- Since A^(p) is the same for all permutations, 
it is only -Fek(p) which distinguishes between them.^° The argument of the logarithm in 
Fek{p) is just the lattice gluon propagator, and the overall factor d — 2 is the number of 
transverse gluons. The key observation is simply that -Fek(p) depends on the permutation 
of the momenta. To see this qualitatively, note that, because the logarithm is a concave 
function, -Fek(p) is minimized by choosing permutations in which there are values of (ab) 
for which Ap'J^ is simultaneously small in all directions. In other words, one lowers the free 
energy by aligning, or "locking", the momenta in different directions. The gain one makes 
by locking the small momenta outweighs the loss incurred by locking large momenta. 

We note that it is the same free energy -Fek(p) that causes the spontaneous breakdown 
of the center symmetry in the EK model. In that case the momenta are fully dynamical, 
and -Fek(p) causes them to be equal, as discussed in Sec. [llj This collapse is prevented by 
quenching, but quenching, which acts independently in each direction, does not prevent 
correlations between momenta in different directions, such as that induced by "locking" . 

We have numerically checked the argument that minimizing Fek{p) leads to locking in 



the following way. We considered the clock momenta, and evaluated F^k of Eq- (4.4) for 



many random permutations of the momenta. What we find is that the vast majority of 

Note that if one uses the weight function pvdM(p) then i^2(p) is canceled by the Vandermonde determinants 
in pvdM- If one uses the clock momenta then A^(p) is a constant. 
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permutations have a free energy larger by 0{N'^) than that for the completely locked case. 
(An example of this result is given below in Fig. Isl) Thus as 6 ^ oo at fixed A^, the 



completely locked vacua dominate. As noted in Sec. HID, the {Z^^Y transformations and 
reflections form a subset of the permutations, and for these the free-energy is invariant. 
Thus there are (2A^)'^~^ degenerate vacua of the locked type, whereas for general (non-clock) 
momenta we expect only a single vacuum. 

In preparation for the numerical study, we now discuss quantities that can be used to 
discern the predicted locking of momenta. As we will explain, for the clock momenta these 
are appropriately called order parameters, although for general p{p) they are not. The 
simplest choices are the expectation values of the d{d — 1) open loops 

M^^, = ti{U^U,)/N and M^,_, = tr (f/^f/t)/Ar (/i > z/) , (4.6) 

which are sensitive to correlations between gauge flelds in different directions. The utility of 
these quantities is particularly clear for the clock momenta, for which one of the permutations 
leads to the A^ being equal in all directions. Then, if — 1, half of the |M^y| equal unity 
(M^ = 1), while the other half vanish (M^ ^, = 0). The same absolute values of the M^ ^, 
hold for the other locked vacua obtained by acting with {Z^Y transformations, while M^_y 
and M^^y switch roles under reflections. This suggests using the combined quantity 

M = 5](|M^,,| + |M^,_,|), (4.7) 

^^<u 

as a signal for locking. We use both M and the individual M^y in our numerical study. 

To illustrate the utility of M, we present in Fig. [3] a scatter plot of the normalized free 
energy versus M for a large set of randomly chosen permutations of the clock momenta 
and with = 1. We include the locked vacua by hand, since they are not among those 
chosen randomly. The flgure indicates that the locked conflgurations have free energies 
that are at least of 0{N'^) smaller than those of the "unlocked" conflgurations, and have 
signiflcantly larger values of M . Taking the results at face value, one might be concerned 
that the number of unlocked vacua might overcome their higher free-energy. However their 
entropy factor is InA^! ~ A^lnA^, which is thus smaller than the free-energy difference of 
0{N'^) or greater. Nevertheless, this plot does lead one to expect that, for flnite A^, a range 
of "partially unlocked" states (present in the flgure only for A^ = 10) will be populated. 

In Section |V] we use these order parameters and other numerical evidence to argue that 
locking does occur also nonperturbatively. For the remainder of this section we discuss in 
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FIG. 3: The dependence of / = F-EK{p)/{{d — 2)N{N — l)/2) on the combined order parameter 
M for random permutations of the clock momenta. Results are for d = 4 and N = 10,20,40,200. 
The fully locked points with M = Q are included by hand. There are also partially unlocked states 
which, for each A^, interpolate between the mass of unlocked states and the locked ones. These do 
not appear in the random sampling except for = 10. 

more detail how locking leads to a failure of reduction. In particular, we explain why the 



previous arguments, summarized in Section III C do not hold. 



B. Implications of momentum locking for large- reduction 

In this section we first focus on the case of the clock momenta, and then return to 
other choices. As noted above, the quenched average in this case is just an average over 
permutations of the momenta. But we now understand that the non-perturbative QEK 
model automatically includes this sum over permutations — it is self- averaging. Thus, in 
principle, the additional quenched average is unnecessary. We also know, however, that 
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the permutations are included with different relative weights — this is manifest in the weak 
coupling free-energy landscape of Fig. [3} and there is no reason to expect equality for other 
couplings. Regardless of the details, the mere fact that the weights are different implies 
that the integrations over momentum space that are induced by the sums over color indices 
are not uniform. This is sufficient to invalidate reduction — the momentum integrations in 
the reduced and infinite- volume cases are different. The case of complete locking provides 
an extreme example: the momentum of each gluon then has the same component, Ap^^, 
in each direction, and the integration over the d-dimensional Brillouin zone collapses to an 
integration along the 1-dimensional body diagonal. 

The argument for reduction based on the loop equations also fails, because one or other 



of the key steps, eqs. (3.11) and (3.12), does not hold. To see how this works, we write out 



these relations for the case that IVopen = M^^^ and PVopen = ^^l,u■ 



I dp (m^,, m;,)^ = Jdp (M^,,)^ (m;,)^ + o(i/iv2) (4.8) 

dp (M^,,)^ (m;,)^ = Jdp {M,,,)^l dp'{M;^,)^, + Oil/N). (4.9) 



We now argue that, if locking occurs, then, for some fi and u, the following two statements 
are correct : 



(I) Ther.h.s. of Eq. (4.9) isof 0(1/A^). 



(II) Thel.h.s. of Eq. (4.8) is of 0(1). 



Thus one or both of the relations must be wrong. The numerical evidence of Section [V] 



suggests that it is the second relation, Eq. (4.9), which fails 



It is easy to see that statement (I) is correct regardless of the choice of p{p). This is due 



to the center symmetry (3.7), under which 



(M,,.)^ (M,,.)^^,^„/^ = (M,,.)^ e2'^^("''+-)/^, (4.10) 

Since the measure dp is unchanged when p — > p + 2iTn/N, the phase factors will cause 
/ dp {M^^J)^ to vanish. For the clock momenta the situation can be slightly different. There, 
self-averaging may take place, and this means that the momenta that contribute to {M^^J)^ 
are all those related to p by permutations. These include also the momenta p^ + 27™^/A^ with 
integer, and so {M^^^)^ ~ I]n^,n„ g27ri(n^+n„) _ Consequently we see that self-averaging 
in the clock momenta case makes statement (I) correct even without the integrations. 
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To see why statement (II) is correct note that the locking means that some of the 



will have 0(1) values. In contrast to the integrands of Eq. (4.9), the integrand 
here, M^^^M*^ = |M^,^p, is invariant under the center symmetry, and thus maintains its 
0(1) value even after integration. 



Which of the two Equations (4.8) and (4.9) fails? This depends on the nature of the 



dynamics. If the self-averaging occurs, then the second step, which for clock momenta is 
simply 

{M,,u), (m;,)^ = (M,,.)^ (m;,)^, , (4.11) 

is trivially valid, and it is the first step which fails. This breakdown of large- factorization 
is then an example of the breakdown of cluster decomposition due the presence of multiple 
vacua — all those related by the center and reflection symmetry. 

The other possibility is spontaneous symmetry breaking (SSB) of the center symmetry, 
in which the system gets stuck in the vicinity of one of the locked vacua. We recall that Z[p) 
itself is {Z^Y symmetric with clock momenta, so there is a symmetry to break. Furthermore, 
despite the fact that the QEK model has zero volume, SSB is possible when N oo because 
there are then an infinite number of degrees of freedom. The (M^ j,)^ (or, indeed, the 
quenched expectation values of any open loops) are order parameters — non-vanishing values 
indicate SSB. If SSB takes place then, by definition, self-averaging no longer occurs, and 
vacuum expectation values of open loops vanish only if we explicitly average over the input 
permutations. If the input permutation is changed by a center transformation, then, since 
the dynamics is {Z^Y invariant, the vacuum that is selected will also be changed by the same 



transformation. In this possibility, factorization, Eq. (4.8), holds, because fluctuations about 
the single vacuum are suppressed as ^ oo. It is the second step, Eq. (4.9), that fails. On 
the l.h.s. the same vacua are selected in the two quenched expectation values, because the 
same input momenta are used, while on the r.h.s. different vacua are, in general, selected. 
Thus the l.h.s. will be of 0(1) for all input p, while the first term on the r.h.s. will average 
to zero. Thus what we call quenched factorization fails. 

We discuss which of the two possibilities — self-averaging or SSB — is expected to occur in 
the next subsection. Regardless of which occurs, however, the key point is that the combi- 



nation of the relations (4.8 4.9) fails, either invalidating cluster decomposition or breaking 



the center symmetry, and thus large-A^ reduction fails. Furthermore, one can numerically 



test for this by calculating the l.h.s. of (4.8) and determining whether it falls as 1/A^ (as 

21 



required for reduction) or tends to a constant as ^ oo (reduction fails). 

We now consider the uniform weight function, p{p) = 1. In this case Z{p) is not center- 
symmetric and the M^^,^ are not order parameters. Nevertheless, if locking occurs, we expect 
something similar to SSB to take place. For a random input choice of p, we expect the system 
to sample the space of permutations, until it finds that with the smallest free energy. Note 
that none of the permutations will be related by center or reflection symmetries, so all are 
expected to have different free energies. In this picture, the system ends up fluctuating in 
the vicinity of a particular permutation. If the weak-coupling free energy is any guide, the 
chosen permutation will be such that if, for a given pair of indices (a, b), the difference Ap'J^ 
is small for one value of /i, then it will also be small for all other values of fi. In other words 
the chosen momenta will be partially locked. Thus, even though the input momenta are 
uniformly distributed in the BZ, those chosen dynamically are not, and planar perturbation 
theory is not correctly reproduced. 

The expected partial locking implies that, for most input p, some of the M^^^, will fluc- 
tuate around complex values with magnitudes of 0(1). If so, this invalidates the quenched 



factorization of Eq. (4.9), because the l.h.s. averages to an 0(1) value, while the p-integrals 
on the r.h.s. implement the center symmetry and cause the averages to vanish. This picture 
is conflrmed by our numerical flndings in Section [V] 

For pbz the situation is similar to that for pciock- There are multiple locked vacua related 
by the center and reflection symmetries, and locking invalidates reduction. The difference 
is that it is only the {Z^Y subgroup of the full center symmetry which is realized, where 
K = N^/'^. To make clear how the presence of permutations in the dynamics unravels the 
carefully chosen coverage of the BZ, we can refer to the simple example in Fig. [2j Permuting 
the momentum components in the "1" direction as, for example, pi ^ p\ for [ah) = (25), 
moves the two momenta in the "boxes" labeled 2 and 5 in the Figure into those labeled 
by 1 and 6, where the momenta are locked. Similarly, all other off-diagonal pairs can be 
moved by permutations onto the diagonal. Thus if the free-energy favors locking, as the 
weak-coupling argument implies, then the momenta chosen by the simulation will lie on the 



The complete locking possible for clock momenta is not possible here because the components of Ap's in 
different directions are different. 
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one-dimensional diagonal of the BZ. 
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Finally, we briefly discuss the choice p{p) = pvdM(p) of Eq. (3.16). This is in some sense 
intermediate between the clock and uniform choices. On the one hand, any value of p is 
possible with pvdM(p), while, on the other, the large- limit of pyduip) is Pciock(p)- Thus 
we expect locking or partial- locking for pvdM, and this is indeed what we flnd numerically. 



C. Expected size of fluctuations 

In this subsection we address the question of whether, for the clock momenta, we expect 
the theory to exhibit SSB or not. We are interested in this question for flxed b and — >■ 



oo. The weak-coupling result of Eq. (4.4) provides a guide to the free-energy landscape, 
and suggests that the dominant states correspond to fluctuations about the locked vacua. 
As in any statistical mechanical system, the issue is whether the fluctuations are large 
enough to cause the theory to move from one locked vacuum to others related by symmetry 
transformations. In inflnite volume, we know from the Mermin- Wagner- Coleman theorem 
[36] that for (i > 2 the fluctuations are not IR divergent and SSB is possible, while for d < 2 
it is not. The question is how this result translates to the QEK model where the spatial 
volume is embedded in the color space. 

To get a rough idea of what happens, imagine that we are in a completely locked vacuum. 
A measure of the fluctuations in the (normalized) traces of open Wilson loops (such as the 
M^^^) is given by the "tadpole" graph 



where p is flxed, and we have used Eq. (3.15 ). The g comes from expanding the t/^, and the 
1/A^ from the normalized trace. Note that since we are doing perturbation theory we can 
really flx the momenta, and we are taking p to be locked. This means that | Ap°^| = \P"' — P''\ 



is independent of z/, and the tadpole can be rewritten as 

N sin2([P« - P^]/2) ' Jc/N 7T ~sm\q/2) 



Tiocked = C E .■.^2nJ mi/o^ = V& r ^—2^, + Oil/N) . (4.13) 



Here we note again that, when any of the p° are equal, as they are for Pbz{p), then there are flat directions 



which are not Gaussian, and the form Fek{p) of Eq. (4.5 1 is invalid. As mentioned above, we do not 
study further the effect of these flat directions, but rather investigate the QEK model with Monte-Carlo 
simulations (see next section). 
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As ^ oo, the sum has gone over to an integral, but the integral is over a one- dimensional 
momentum space, and is thus IR divergent. The cut-off c/N (with c a constant that could be 



determined by a more complete analysis) arises from fact that the original sum, Eq. (4.12), 
has a minimum Ap of 0(1/A^). The IR divergence implies that 

TiockedCxAr/6[l + 0(l/Ar)] (4.14) 

so that the fluctuations about the locked vacua diverge as ^ oo for fixed b. 

These divergences can be anticipated from the result for the maximum energy barrier 



^'S'qek (see Eqs. (4.1)-(4.2)) that exists between two configurations related to each other 
by the permutation ^ P^i- Denoting Apl^ = pIj, ~ "we see that if Ap^^ = 0{1/N) and 
Ap'^ ~ 0(1), then IASqekI ~ b/N , and fluctuations in the direction parameterized by the 
SU{N) matrix (4.1) overcome the barrier when A^ is large enough that b/N < 0(1). For 
locked vacua with Ap^ = 0{1/N) in all directions, the barrier is even lower, | ASqekI ~ b/N^. 
The Gaussian terms in the action for these "flat" directions is small relative to higher-order 
terms, and ignoring the latter in the tadpole calculation leads to the apparent IR problem. 

It is, in fact, the more severe b/N^ divergence which leads to Tiocked oc A^. One can see 
this by noting that for a random permutation of the clock momenta T ~ / df'q/q^ in the 
IR, and this is convergent. This is despite the fact that there are the flat directions with 
|A5qek| ~ b/N. 

The upshot of this discussion is that we cannot quantitatively trust the weak-coupling 
calculation of T for the locked vacua if A^ ^ cxd at flxed b. How does this affect the free- 



i^EK(piocked) ~ Ar2 ' dqq\og{q) ~ 0{N^) + 0(Ariog AT) , (4.15) 

Jc/N 



energy Fek{p) which we discussed above for the locked vacua? It follows from Eq. (4.4) 
that 

rO{l) 
I c/N 

and so the leading order term is IR safe, while the subleading term cannot be trusted in a 
Gaussian analysis. Since our previous discussion was based on the leading 0(A^^) term, it 
remains valid. 

Returning to the issue of SSB, we need to know whether the large- A^ divergence of Tiocked 
implies that the system will fluctuate into nearby locked vacua (which have momenta differ- 
ing by center-transformations or reflections). We know that there will be large fluctuations 
in the directions given by transpositions between close momenta, for these are the source 
of the IR divergence. Thus to address the question we proceed as follows. It is possible to 
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move from one locked vacuum to another by stringing together a sequence of transpositions 
involving nearby p's (i.e. with Ap always of 0{1/N)). As we proceed along such a string, 
the momenta become partially unlocked, and the energy barriers to transpositions increase 
from of 0{1/N^) to of 0{1/N). Nevertheless, they all still vanish as ^ oo, so there is 
a vanishing energy barrier between locked vacua. What matters, however, is whether there 
is a free-energy barrier. We can investigate this using the weak-coupling result, evaluating 



Eq. (4.4) numerically for each momenta along the path. 
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FIG. 4: The free energy F^]^(ji)/[d — 2), divided by N (left panel) or A^^ (right panel), along the 
two paths between locked vacua described in the text. Results are for d = 4 and N = 10, 50, 
200 and 1000. The horizontal axis gives the fraction of the total transpositions required, with 
the starting point being plotted at position 1/A^ (left panel) and 2/[N{N — 1) + 1] (right panel). 
The left panel is for path 1 (vacua related by center-symmetry) , the right panel for path 2 (vacua 
related by a reflection). 



We have considered two types of path. Both begin from a locked state in which t/^ = Ai 
for all /i, where (Ai)^^ = exp(iP"), with the P'^ given in Eq. (3.17). Thus, for example, 
^^,1/ = and Mn^-j^ = 1. Path 1 arrives at a state with Ui = U3 = f/4 unchanged and 
U2 = so that M12 = 0, Mi,_2 = e+^^^^Z^, etc.. The path is made of a series of 

A^ — 1 transpositions between adjacent indices, for each of which | Ap| takes its minimal value 
of 27r/N. The paths we use are exemplified by the following sequence for A^ = 6 (which 
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shows the ordering of the momenta P°, a G [1,6], along the diagonal of U2) 



123456 213465 312564 412356 512346 ^ 612345 . (4.16) 

In contrast, path 2 takes us to a vacuum with Ui = = Uj^ = Ai and U2 = A|, for 
which Mi^2 = 1 and Mi^„2 = 0. This can be achieved with transpositions alone along a more 
complicated path of length A^(A^ — 1)/2, that for = 6 would be the string of transpositions 



in Eq. (4.16), followed by 



612345 6213465 631245 ^ 641235 651234 652134 ^ 653124 ^ 654123 
654213 ^ 654312 654321. (4.17) 

We show the results in Fig. |4j We find that, for path 1, the free energy barrier AF scales 
asymptotically with A^, while for path 2 it scales with A^^. Since we know from above that 
the 0{N'^) part of the free energy is IR safe, we conclude that fluctuations along path 2 are 
certainly suppressed. For path 1 the situation is more subtle, as we now discuss. 

The issue for path 1 is whether the leading contribution to AF, which we see to be of 
0{N), is IR safe, given that F itself is untrustworthy at this order. The numerical results 
themselves suggest that the A^logA^ terms cancel in AF, but it would require a more 
detailed analytic analysis to demonstrate that this cancellation of untrustworthy terms is 
itself trustworthy. Thus the most conservative conclusion is that we do not know whether 
the barrier path 1 grows with A^ and suppresses fluctuations. Other uncertainties in this 
analysis are that we have not investigated all paths, nor accounted for a possible entropy 
factor involving the number of paths, and finally that it is based on the leading term in the 
weak-coupling analysis. Thus to learn about the extent of locking, and the possibility of 
SSB, we must study the QEK model non-perturbatively, and to this we now turn. 



V. NON-PERTURBATIVE LATTICE STUDY 



In this section we present our numerical results for the QEK model. In Sec. VA we 



briefly describe the methodology, focusing on an explanation of the two strategies we adopt 



to perform the quenched average: self- averaging and explicit averaging. In Sec. |VB| we map 
the phase structure as a function of the bare coupling, b, using measurements of the plaquette 
and the order parameters M^^j^. These results lead us to investigate various features of the 
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model in more detail. In Sec. VC we describe results from high-precision measurements of 
the plaquette. This allows us to study the dependence of the results on the calculational 



strategy and on the choice of the weight function p{p) (defined in Section IlIC). We present 



similar measurements for M^^^ in Sec. VD and use them to understand the structure of the 



vacua of the QEK model. Finally, in Sec. |V E[ we analyze a "strong-to-weak" transition that 
occurs in the model, using an adaptation of the Wang-Landau algorithm [37] to perform a 
precise measurement of the coupling ht at which it occurs. 



A. Methodology 



The QEK model has been defined above in eqs. ( |3.1 )-(3.4). The ingredients for a sim- 
ulation are a weight function for the momenta, p(p), and the coupling h in the quenched 



action, eq. (3.2). Specifically one is instructed to draw momenta weighted by p(p), construct 



the diagonal eigenvalue matrices using eq. (2.13), and then do a Monte-Carlo average 



over the SU{N) matrices for fixed A^. Observables involving gauge links, such as the 
plaquette, can then be reconstructed using the definition = V^A^V^. 

As noted above, we consider four choices of weight function: uniform (puniform(p) = 1); 



clock [defined in eq. (3.20)], Vandermonde [defined in eq. (3.16)], and BZ [defined above 
eq. ( |3.21 )]. It is straightforward to draw momenta from the first three of these distributions, 
while there is only a single choice for the BZ distribution. 



The Monte-Carlo integration over the is non-standard because the action (3.2) is 



quartic in each of these matrices, so that a simple heat-bath algorithm cannot be used. 
Instead, we use the following three approaches. 

1. A straightforward (and slow) Metropolis algorithm using the original action, updating 
all of the SU (2) subgroups of each in turn; 

2. A faster Metropolis algorithm, using a Gaussian auxiliary field to reduce the action to 
quadratic order in the [38]. We again update SU{2) subgroups of in turn. 

3. A combination of a Cabibbo-Marinari heat-bath (again applied in turn to SU{2) sub- 
groups) and various type of over-relaxations (both SU{2) and SU (N), the latter using 
the method of Ref. [39]). These are applied after using two Gaussian auxiliary fields 
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to make the action linear in the V^. We use a ratio of one heat-bath update for every 
four over-relaxations. The details of this algorithm are described in Ref. 



We find that the second algorithm typically decorrelates our measured quantities most 
rapidly, and we use this for most of our runs. All the measurements in this paper were 
separated by 5 full updates of all four V^'s. 



We perform the evaluation of the quenched average, Eq. (3.4), using one of the following 



two strategies. The only exception is for pez, for which no average is necessary. 
Strategy A : Explicit quenched averaging 

Here we simply follow the quenching recipe laid out above: generate an ensemble of sets 
of momenta weighted by the distribution p(p), calculate {0(U))p for each member of this 
ensemble, and then average over the ensemble. We have analyzed the QEK with this strategy 
for all choices of p{p) listed above, except pbz, but have mostly focused on Pciock(p)- 

Strategy B: Self-averaging 



As mentioned in Sec. |III C| if reduction is valid, then one need not perform the momentum 
integral at large- — a single value p = po is sufficient, since the sum over color indices will 
sample the Brillouin zone. We refer to this possibility as self-averaging. For the choice 
p{p) = Pciock(p), self-averaging is, in principle, exact for all finite A^, because, as explained 
in Sec. [nTD| 

(0(W))qek = (0(W))po • (5.1) 



Here, for each p, (po)/^ can be any permutation of the clock momenta Pa defined in Eq. (3.17). 



We recall that Eq. (5.1) holds because the integral over includes all permutations of the 
elements of A^. For p = pciock, we often use {p'^)o = -P" for all p (which we call po = Piockod)- 
The self-averaging strategy is not guaranteed to work in practice, because it may be 
that the simulation fails to fully explore all possible permutations due to algorithmic short- 
comings [221 SO]- To check whether this happens it is useful to measure quantities which 
change as one moves from one permutation to another, and we use the order parameters 



Mf^^^ [Eq. (4.6)] for this purpose. 

It is important, however, to distinguish such an algorithmic failure from a genuine break- 
down of reduction. In the latter case, most of the permutations will have higher free energy, 
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and will be visited with a probability which vanishes as ^ oo. Furthermore, if sponta- 
neous symmetry breaking occurs (as is possible with Pciock or pez) then, as ^ oo, the 
system fluctuates around a single vacuum due to the infinite barrier between vacua connected 
by particular permutations. 



B. Mapping b dependence 

We begin by determining the dependence of the average plaquette, 

1 



V 



ReJ2{TTU,UMlUl), (5.2) 



Nd{d -l)/2 

and of the M^^^,, as a function of b in the range 0.1 ^ h < 1. We use the self-averaging strategy, 
and have used all choices of p{p). Our focus in this section is on qualitative features, and at 
this level we do not find much dependence on the choice of p{p). For brevity, therefore, we 
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only present results for p{p) = Pciock- 

We use "hysteresis runs" starting either from a "cold" field configuration with 

V^ = l, (5.3) 

at a high value of b = 6coid — 0.7 — 1.0, or from a "hot" field configuration with 

= random element of SU (N) , (5.4) 

at a low value of 6 = 6hot — 0.1 — 0.2. For a cold (hot) start we gradually decrease (increase) 
b until we reach the value b = 6hot(coid)- We study gauge groups with 10 < A^ < 200, and list 
the parameters of our major runs in Table [!} 

In Fig. [5] we present results for Up from simulations with A^ = 50, 80 and 100. At first 
glance, the plots appear consistent with the validity of reduction. The results for Up are 
close to the analytic predictions and to the numerical results from large volume, large- 
A^ simulations (although we do not show the latter here). The increasing hysteresis with 
increasing A^ is indicative of a strongly first order phase transition somewhere in the range 
bt ~ 0.30 — 0.35, which is indeed close to the couphng, 6buik — 0.36, at which the well-known 
"bulk" transition occurs in the large- A^ gauge theory 



We do see, however, that the 1/iV corrections for /Ouniform(p) are very large for b < 0.3 and smear a 
strong-to-weak transition that occurs at about b ~ 0.32. This observation was also reported in Ref. [T5]. 
There is also an estimate of 6buik from the TEK model [42], the equivalence of which to large-TV pure 
gauge theory has been thrown into doubt by the work of Refs. [371 ISHl ■ 
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N 



equilibration updates measurements 



p(p) and p 



20, 40 
20, 40 
50 

80 

100 

16 

81 



100 
1000 
1000 

100 

500 
1000 
100 



200 
5000 
5000 

200 

1000 
5000 
500 



uniform, VdM 

clock (piocked) 
clock (piockod) 

uniform 

VdM 

clock (blocked) 
clock (piocked) 

BZ 
BZ 



TABLE I: Details of hysteresis runs used to map the phase diagram. For each value of b, we use 
the quoted number of equilibration updates followed by the quoted number of measurements (the 
latter being made every 5 updates). 

Below, and in the next three sub-sections, we show that this impression is wrong. A 
clear signal for this can be seen in Fig. ([6]), which shows how the absolute values |M^ ,^| 
depend on b. We recall that the transform non-trivially under center and reflection 



symmetries. The discussion of Sec. IVB implies that, for reduction to hold, the |M^^j,| 



must have expectation values of 0{1/N), and thus that the (|M^^j,|) should fall to zero as 
— > oo. In fact what we find is that some of the M^^^^ fluctuate around 0(1) values in the 
weak-coupling phase. This is the first indication that reduction does not hold in the QEK 
model. 

This result calls for a more detailed study of systematic errors. These include the pos- 
sibility of very large 0(1/A^) corrections (i.e. that the non-zero values for (|M^^^|) would 
decrease for large enough A^), dependence on the choice p{p) or on the self-averaging strat- 
egy, and the possibility that the simulations did not, in fact, equilibrate (and that given 
enough updates would tunnel into a "vacuum" that satisfies reduction). In addition, a more 
accurate determination of the transition coupling bt would allow a direct test of reduction. 
This is the coupling at which the QEK model goes through a first order transition, and, if 
reduction holds, should equal the bulk-transition coupling febuik of the infinite- volume, large- 
A^, pure-gauge theory. In the past, the numerical proximity of bt and 6buik was considered 
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FIG. 5: Hysteresis plots of the plaquette variable Up versus b for SU{50) ([blue] crosses), SU{80) 
([magenta] open squares), and SU{WO) ([light blue] filled squares). Results are for Pciock and self- 
averaging. The curves are the predictions for SU{oo) in the strong-coupling expansion to leading 
order (solid [red] curve) and of the weak-coupling expansion to three loop order (dashed [blue] 
curve) (taken from, for example, Ref. [26]). The lower-panel shows a close-up of the strong-to- 
weak transition region. 

31 



















hot Stan 
cold start 


















' 1 ' 


1 1 = 


^ 1 -! 
































t ^ 






; : 




















\ \ 




















1 1 


































i ' 












1 « • 


1 * 












! i i 


1 1 1 


i i , 



0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 



















hot Stan 
cold Start 


1 - 
















« * 

» 


« 1 

( 










1 


I 

i 

1 








































i i 


( 


























































■ * * 


< •< 


! 


i 


( i 




» • « 









0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 



FIG. 6: The values of all twelve (|M^jy|) plotted versus b for A'" = 20 (upper panel) and A'" 
(lower panel), using pciock and self-averaging. 
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choices of p equilibration updates measurements p{p) 



20, 40, 80 
40 



20 
20 



1000 
5000 



5000 
1000 and 2000 



clock 
uniform 



TABLE II: Simulation details for measurements of Up using strategy A — explicit quenched averag- 
ing. The number of equilibration updates and measurements is for each value of p. We use b = 0.4, 
0.45 and 0.5 in all cases. 



N 



equilibration updates measurements 



choice of p{p) and p 



20, 40, 80 

20, 40 
80 
50 

100,125,150,200 
16, 81 



5000 

5000 
1000 - 5000 
1000 
500 
5000 



100000 

100000 
10000 
20000 
1000 
100000 



clock {p = Plockcd 

or permutation) 
uniform 

uniform 

clock {p = piockod) 
clock {p = Plockcd) 

BZ 



TABLE III: Parameters of simulations used to calculate Up using strategy B — self-averaging. We 
use b = 0.4, 0.45 and 0.5 except for N = 50, where we only use b = 0.4 and 0.5, and for = 100, 
125, 150, 200, where we only use b = 0.4. 

to be evidence in favor of large- equivalence [HI [121 EH [13 HB] ; but the calculations of bt 
were not of high accuracy. In the next sub-sections we attempt to address all these issues. 



C. Precise measurements of the plaquette 

In this section we perform high-precision measurements of Up for b = 0.4, 0.45 and 
0.5, values chosen to allow comparison with results from the large volume simulations of 
Ref. [271111]. We use and test both strategy A (explicit quenched averaging), and strategy B 
(self-averaging), and in addition study different choices for the measure p{p). To implement 
strategy A we draw a new choice for p (drawn randomly with weighting p{p)) after a fixed 
number of equilibration and measurement sweeps. The simulation parameters are given in 
Table [11} For strategy B we simply use very long runs with a fixed choice of p. Details are 



given in Table III 
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We begin by comparing the plaquette time histories for the two strategies. In Fig. [7] we 
show results for = 40 at 6 = 0.40 and for = 80 at 6 = 0.50, in both cases using Pciock(p)- 





» +t * + >4.*i + -f*?^ 



; +f "iln- VF : V ^ +4- 



10000 20000 30000 40000 50000 60000 70000 80000 90000 10000C 



10000 20000 30000 40000 50000 60000 70000 SOOOO BOOOO 100000 





1 1 


1 1 


1 1 


1 1 1 1 

(1.2) 

+ 

+ 


T 

^laquette • 














10000 20000 30000 40000 50000 60000 70000 80000 90000 10000C 




10000 20000 30000 40000 50000 60000 70000 SOOOO 90000 100000 



FIG. 7: The (fi = 1, v = 2) plaquette variable vs measurement number. Equilibration updates 
are not shown. The upper two panels are for 5'C/(40) with b = 0.5, the lower two for SU{80) 
with b = 0.4. The left-hand plots are for strategy A, in which a random permutation of the clock 
momenta is generated every 5000 measurements. The right-hand plots are for strategy B, with 

P = blocked- 

The results for = 40 suggest that, once equilibrated, both strategies give similar 
results. The "tails" below the main band for strategy A indicate, however, that insufficient 
equilibration sweeps were included. This effect is much clearer for = 80, because of the 
smaller fluctuations. To see what happens with sufficient equilibration, we present in Fig. M 
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the time histories for 5'f/(40) at 6 = 0.50, obtained using strategy B, with both p = piocked 
and a single fixed random permutation. After a long period in a metastable state, with a 
relatively low value of Up, the system does appear to tunnel into a state with a plaquette 
value consistent with that for p = piocked- The time this requires (about 15000 measurements) 
exceeds, however, that used in our strategy A runs. 



(1 ,2) Plaquette starting from a locked vacuum 
(1 ,2) Plaquette starting from a random permutation 




0.68 



0.69 



Q 5000 10000 15000 20000 25000 30000 

FIG. 8: The = (1,2) plaquette variable for SU{40) at b = 0.50 vs. update with strategy B 

(self- aver aging) with Pciock- [Red] plusses are for p = piocked; while [green] crosses are for a random 
permutation of the clock momenta. 



We tentatively conclude that self-averaging works at least approximately for the plaque- 
tte, given long enough equilibration times. This conclusion is supported by the results from 



the other simulations listed in Tables O and III To illustrate this, we collect, in Tables IV 



and|V) the average values of Up using both strategies and for clock and uniform densities. For 
strategy B we show only results from runs that were equilibrated. For strategy A all results 
are suspect because of the thermalization issues discussed above and illustrated by Fig. [7} 
We nevertheless include them as a comparison. The table also includes the best estimates for 
the plaquette values for SU (oo), obtained by extrapolating from large- volume simulations 
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b, Up for SU (oo) 


Strategy and p{p) 


SU{20) 


SU{40) 


SU{80) 




A, clock* 


0.7294(2) 


0.7256(2) 


0.7223(4) 




B, clock (p = jDlorkpd) 


0.7396(5) 


0.7425(2) 


0.7429(1) 


b = 0.50, Up ~ 0.7182 


B, clock {p = random) 


0.739(1) 


0.7424(1) 


0.7429(1) 




A, uniform* 


0.7401(3) 








B, uniform 


0.7483(3) 


0.7405(2) 


0.7432(1) 




A, clock* 


0.6968(5) 


0.7090(1) 


0.6867(8) 




B, clock (p = Plnrkpn) 


0.7035(9) 


0.7094(1) 


0.7100(1) 


b = 0.45, Up ~ 0.6795 


B, clock (p = random) 


0.703(1) 


0.7095(1) 


0.7101(1) 




A, uniform* 


0.7082(4) 








B, uniform 


0.7153(2) 


0.7086(2) 


0.7100(1) 




A, clock* 


0.6533(5) 


0.6489(7) 


0.645(1) 




B, clock {p = blocked) 


0.66014(54) 


0.6651(2) 


0.6665(1) 


b = 0.40, Up ~ 0.6259 


B, clock {p = random) 


0.6595(5) 


0.6645(3) 


0.6665(1) 




A, uniform* 


0.6648(5) 








B, uniform 


0.6737(5) 


0.6652(2) 


0.6642(3) 



TABLE IV: Comparison of plaquette expectation values between averaging strategies and different 
choices of p{p). The results from strategy A are denoted by a "*" to indicate that they are suspect 
due to a possible lack of equilibration (see text). The first column includes the estimates for SU (oo) 
based on extrapolations using large-volume simulations from Ref. 



j. These are the numbers that would be reproduced by a large extrapolation of QEK 
results were reduction to hold. 



We first comment on the results using strategy B. We first note (from Table IV) that, in 
all cases, the final plaquette is independent of the choice of input momenta for p = pciock- 
This is the expected self-averaging for a center-invariant quantity. More striking is that, as 

increases, results from strategy B using Puniform appear to converge to those from Pciock- 
This gives us confidence that we are not observing systematic errors due to the choice of 
p, and that the systematic differences with the results from strategy A are due to lack of 
equilibration of the latter. 

The most important comparison is with the results for the infinite-volume SU{oo) theory. 
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SU{50) 


SU{100) 


Sf/(125) 


SU{150) 


5[/(200) 


0.6662(9) 


0.6647(3) 


0.6658(3) 


0.6667(3) 


0.6670(2) 



TABLE V: Additional results for Up obtained with strategy B at 6 = 0.40 with /Ociock and p = piocked- 

To make this more precise , we extended the results at b = 0.4 up to = 200 (see Table [V|) . 
The resulting comparisons are shown in Fig. ([9]). We have plotted Up versus 1/A^, since 
this is the expected A^ dependence in the QEK model. Our results show a fairly smooth 
extrapolation to A^ = oo, with small corrections whose dependence on 1/A^ we cannot 
definitely determine. We do not perform a detailed fit, however, since it is clear that, 
regardless of the precise form of the subleading terms, our results extrapolate to significantly 
higher values of Up than those of the infinite-volume lattice gauge theory. This discrepancy 
clearly shows that the QEK model does not reproduce the physics of the large-A^ gauge 
theory. 

P^Pclock 

0.76 

0.74 

0.72 

B 0.7 
cr 

E 0.68 
0.66 
0.64 
0.62 

-0.01 0.01 0.02 0.03 0.04 0.05 0.06 

1/N 




FIG. 9: Up versus 1/N, for b = 0.4, 0.45 and 0.5, compared to the expectations for infinite-volume 
SU (oo) gauge theory (presented inside the ellipse). Lines are only to guide the eye. Results shown 
are for p = pclock with p = piocked- 



We have also obtained results using pez- These are collected in Table VF and the compar 



ison to the lattice large- A^ result is shown in Fig. 10 The discrepancy with infinite- volume 
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N 


6 = 0.4 


6 = 0.45 


6 = 0.5 


16 


0.88627(5) 


0.89935(4) 


0.90961(3) 


81 


0.812323(2) 


0.83546(1) 


0.85291(1) 



TABLE VI: Results for Up using strategy B for pBz- 



SU (oo) values is significantly larger in this point we return to below. 
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FIG. 10: As in Fig. |9] but for pez- Note that the vertical scale differs from Figj9j 



D. Precise measurements of the ^ 

To elucidate the nature of the breakdown of reduction, we present here results for the 
"order parameters" M^^^. We use the same simulation parameters as in the previous section. 
We recall that, for reduction to hold, (|M^_,^p)p should be no larger than 0(1/A^) for all \x 
and V. Furthermore, for the case of Pciock, the expectation values {M^i,)p are true order 
parameters for spontaneous breakdown of the center symmetry. 



We begin by presenting, in Fig. 11, the Monte-Carlo time history of the real parts of a 
selection of the M^^^,, using p{p) = pciock(p) and strategy A (explicit quenched averaging), for 

= 40 and b = 0.50. This is the run for which we have previously shown the plaquette in 
the upper-left panel of Fig. [7| We clearly see equilibration into distinct "vacua" for different 
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10000 20000 30000 40000 EOOOO 80000 70000 80000 90000 10000C 



FIG. 11: Real part of M^^^ versus measurement number for = 40 and h = 0.50 with Pc\ock{p)- 
The figure shows a sequence of 20 Monte-Carlo runs, each with 1000 equilibration updates (not 
shown) followed by 5000 measurements, and each with a randomly chosen permutation of the clock 
momenta. For clarity, we only present M2,i ([red] plusses), M4^3 ([green] crosses), and M2,-2 ([blue] 
fancy crosses). 



choices of input momenta, and in several cases we can see the tail-end of what appears to 
be a tunneling process. The values of Re(M^ ,^) either oscillate around zero or around values 
of 0(1). The latter indicate SSB of the center symmetry, and the presence of the locked 
momenta discussed in Sec. IIV A[ 



We show a similar plot for Puniform(p) in Fig. 12, except that we use only a single random 
choice of p, have longer runs to assure equilibration, and present results for both = 40 
and 80. After a long equilibration period, the M^^y at both A^ fluctuate around what we 
assume to be vacuum values. We note that the fluctuations are smaller for the larger A^, as 
expected in general. We see this behavior throughout our study. The crucial observations, 
however, are that some of the M^^j, fluctuate around non-zero 0(1) values (indicating locked 
momenta), and that these values are comparable for both A^. This implies that the left-hand 



side of Eq. (4.8) is of 0(1), and reduction does not hold. 

It is also instructive to look at the full complex values of the M^^^. In the left panel 



of Fig. 13 we show the scatter plot for the same data-set used in Fig. 11, Apart from 
"equilibration tails" , we see that the simulations settle down into vacua in which a given 
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SU(40), M(a2), M(4,1), M(4,2) 




-0.6 I ' 1 ' 1 1 

2000 4000 6000 8000 10000 



FIG. 12: Real part of M^, ^ versus measurement number for Puniform 

(p), with a single random p, 

at 6 = 0.5 and for = 40 ([red] plusses) and 80 ([green] crosses). We show only M3^2, M^^i, and 
M4,2 for SU {40), and M4,2, M4,3, and Ah,-! for SU (80). 




plot of all twelve M^^jy obtained from runs with p{p) = /Juniform(p) for SU{40) ([red] plusses) and 
SU{80) ([green] crosses) at 6 = 0.50, with only equilibrated results shown. 



either fluctuates around or around mo exp(27rm/40), with n an integer and tuq ^ 0.65. The 
different vacua are related by (an appropriate subgroup of the (Zjv)^) transformations. This 
is quahtatively consistent with what we would expect with locked vacua when fluctuations 
are included. Without fluctuations, the locked vacua have half of the M^,^ vanishing, and 
the other half of the form exp{27cin/N). The fluctuations reduce the magnitude from unity 
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to mo- Note that this reduction is greater than one would predict from a simple mean link 
model, in which rriQ ^ ^ 0.86. This may be a consequence of the fact that the partial 
unlocking of momenta can reduce \M^^\ while leaving the plaquette unchanged. 

This figure gives a very clear illustration of the way in which (M^^,y)QEK vanishes when 
using pciock(p)- One is instructed to average over input momenta which are permutations 
of the clock momenta. For given input momenta, the dynamics picks a (partially) locked 
vacuum. As the average is taken, each will end up with equal probability in the center 
near the origin, or in the "ring" of radius tuq, and in the latter case with equal probability 



in each of the N vacua. In this way M^ ,^ will average to zero. As noted in Sec. IV, the 
dynamics will determine whether, for a given input momenta and as ^ oo, the theory 
gets trapped in a single vacuum or moves between them. Our numerical results strongly 
indicate the former, in which case (for pciock) SSB is occurring. 



A similar scatter plot for Puniform is shown in the right panel of Fig. 13 In this case all 
twelve M^ t, are shown for each A^ (not just the three for each A^ shown in Fig. 12), and 
we display only measurements after equilibration. For Puniform there is no center symmetry, 
but we do see (most clearly for A^ = 80) the expected pattern for locked momenta of six 
non-zero and six near-zero magnitudes. (Note that some of the [red] A^ = 40 points near 
the origin are obscured by the [green] A^ = 80 points.) We also observe no reduction in the 
0(1) magnitudes as A^ increases from 40 to 80 — indeed the magnitudes seem to increase. 
This we take as strong evidence for the breakdown of reduction. 

Finally, we consider pbz- We show results obtained only from a hot start. In the 



left panel of Fig. 14, we show the time history of all the M^ ,^ for A^ = 16 and b = 0.40. 



Recalling the definition of pez from Eq. (3.21), we note that, since K = vl6 = 2, all 



are either or tt. This means that the M^^^, are real, and that M^^^ = M^^^^, (so there are 
only 6 independent M^^^,). Furthermore, the center symmetry is only {Z2Y, although this 
symmetry group is still sufficient to forbid expectation values for the M^^^. What we see 
from the figure is that while four of the fiuctuate around zero, two of them (Mi 2 and 
M4^2 ) acquire nonzero expectation values that break the (^'2)^ symmetry. 

We can understand this pattern of expectation values in the following way. The input 



The fluctuations in the runs beginnings from cold starts were too small to allow the simulation to forgets 
its initial state, be it a state with zero or nonzero 



41 




10000 20000 30000 40000 50000 60000 70000 80000 90000 10000C 



0.07 
0.06 
0.05 
0.04 
03 
02 
01 


-0 01 

-0 02 
-0 03 
-0.04 



1 1 


1 1 1 1 

... .A 


1 


T 










i ' " 


: :Xr- 


















w 






i- - - 














- A 








1 - - i - - ■■■■ 




♦ 


- - L ^- - 4 - - ■- 







FIG. 14: Results for M^^j^ with pBz- Left panel: Time history (versus measurement number) of all 
six independent M^^i^ for = 16 and b = 0.40. Right panel: Scatter plot of all Mfj_^,y for = 81 
and b = 0.70. 



momenta [defined in Eq. (3.21)] are such that 



Ai = diag(cr3,o-3,. . .) , A2 = diag(l2, -I2, 12, -I2, • ■ 
A3 = diag(l4, -14,14,-14), A4 = diag(l8, -Ig) , 



(5.5) 



where 1„ indicates (the diagonal part of) an n-dimensional identity matrix. With these 
matrices, and assuming f/^ ^ A^ (i.e. ignoring fluctuations due to the V^), all the M^^j, 
vanish. By a single transposition, however, one can change A2 to 



A2 = diag(-(j3, -I2, 12, -I2, 12, -CTs, I2, -I2 



(5.6) 



One then finds, in the same approximation of ignoring fluctuations, that there are two non- 
zero Mo^j^: Mi^2 = ^^2,4 = —0.25. Fluctuations will reduce the average from this value. Thus 



this scenario provides a possible explanation for the results of the left panel of Fig. 14 This 
is, in fact, one of many choices of transpositions that leads to this pattern of expectation 
values. Furthermore, all the patterns of values for the M^^^, that we have observed in our 
= 16 runs can be explained similarly. 

We see an analogous phenomenon for = 81. Here, since K = 3, the center symmetry 



is (Zs)^. The right panel of Fig. 14 shows a scatter plot of all the (now twelve) M^^u from a 
simulation at 6 = 0.70. One can understand this figure by calculating the possible values of 
M^^i, that are obtained by permuting the elements of the initial A^, and ignoring fluctuations. 



-0 06 -0 04 -0 02 02 04 06 
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The result is shown in Fig. 15, and is clearly a good description of what we see in the right 
panel of Fig. 14 
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FIG. 15: Some of the possible values for M^^jy in the = 81 case, assuming no fluctuations, i.e. 
when the are permutations of pbz- Note that, in a given simulation, one expects only some 
(at most 12) of these values to be realized. This figure should be compared to the right panel of 
Fig.lUl 



We note that, unlike for pciock, the BZ weight function does not lead to complete or nearly- 
complete momentum locking. In a completely locked state, all the are equal up to center 
and reflection transformations, and this leads, in the example of = 2, to all six independent 
M^j^^y being close to ±1. To reach such a locked state requires many transpositions, however, 
and our results suggest that only a few transpositions have occurred. 

In summary, the numerical results presented in this sub-section indicate that some of the 



"order parameters" M^y acquire 0(1) expectation values, which, as described in Sec. IV 



is inconsistent with large- A^ reduction for the QEK model. For the weight functions pciock 
and Pbz, the expectation values for the spontaneously break the center and reflection 
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symmetries. This breakdown is not apparent in the simplest open loops, i.e. (trf/^) with 
n < N, but is exhibited by more complicated objects like the "corner" variables For the 

uniform and clock distributions, the actual values of the M^^^ are qualitatively consistent 
with the "momentum-locking" predicted by the weak-coupling analysis. That analysis, 
however, could not determine whether the symmetry-breaking or cluster-decomposition- 
violating scenario would hold. Our numerical results clearly favor the former. 

E. Precise measurements of the transition coupling bt 

The plaquette data in Fig. ([s]) strongly suggest that the QEK model has a first order 
phase transition for b somewhere in the range 0.30 — 0.33. This was already noted in the 
early QEK literature [HI [T21 [131 E], and the transition was assumed to be the same as 
that which occurs in the SU (oo) gauge theory at ^Buik — 0.36 (the "bulk transition"). The 
~ 10% discrepancy was attributed to 0{1/N) corrections and/or other systematic errors. In 
this section we revisit this issue, and, in particular, attempt to greatly reduce the systematic 
errors in the determination of bt. 

The main source of uncertainty is the strongly first-order nature of the transition, and the 
consequent metastability. The strength of the transition is indicated by the size of the jump 
in the plaquette, which is ~ 0.3. Although, strictly speaking, there is no transition unless 
oo, already for = 50 there is a significant hysteresis regime of width Ab ~ 0.05, 
and this width increases with A^. Thus an estimate of bt from Fig. |5]has an 0(15%) error 
at A^ = 50, and this error too increases with A^. It does not help to calculate Up on a denser 
grid, because of the metastability. 

One way forward is to use re-weighting, making use of those values of b where tunneling 
between phases occurs. We expect the tunneling probability to fall exponentially with A^^ 
(which counts the number of degrees of freedom and thus is like the volume), and our 
results are qualitatively consistent with this. We find that we can successfully use standard 
Ferrenberg-Swendsen (FS) re-weighting [H] for A^ = 20 and 30, but for A^ larger than about 
40 the method fails because tunneling ceases. 

To proceed we need a method which encourages tunneling. We chose to use the "Wang- 

The breaking pattern depends on the extent of locking. For complete locking, and Pciock, the breaking is 
Zfj — !■ Zn, where the remaining symmetry is the diagonal Z^. 
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Type of re-weighting 


SU{20) 


SU{30) 


SU{40) 


SU{50) 


Ferrenberg-Swendsen 


0.29598(5) 


0.30545(5) 






Wang-Landau 


0.29544(37) 


0.30569(17) 


0.30968(20) 


0.31121(19) 



TABLE VIL Values of the strong-weak transition coupling bt, obtained for /9ciock(p) withp = piocked, 
with two different re- weighting methods. 

Landau" re- weighting method, developed recently in the field of statistical mechanics [TF] . 
This required adapting the method from spin-systems to gauge theories, as well as developing 
a systematic way of estimating errors. Presenting this analysis is beyond the scope of 
this paper, and is presented in Ref. |40j. We note only that this is an adaptive method 
of determining the density of states, which includes a feature that forces motion through 
configuration space. 

We have carried out these calculations only for p = Pciock and with input momenta being 
locked. Since the algorithms are designed to ensure ergodicity, however, we expect that the 
simulations will explore multiple permutations of the momenta, i.e. will be self-averaging. 
Evidence in support of this expectation is that we do see many tunnelings between the weak 



and strong phases for all A^. Table VII gives our results for the transition coupling (defined 
as the peak in the susceptibility). We find that the results from both techniques agree (when 
both are available) despite the very small (0.02% — 0.1%) statistical errors. 



We plot these results versus 1/A^ in Fig. (16). For comparison, we include estimates of bt 



from the old numerical studies in Refs. [HI [T2l [131 E]) as well as the most recent estimates 
of the coupling 6buik at which the bulk transition occurs in the infinite volume SU{oo) gauge 
theory [H]. While our new results are consistent with the old numerical studies of the 
QEK model, it is very unlikely that they extrapolate to the vicinity of 6buik- We stress, 
however, that to make this observation, it is crucial to have very small errors, and this was 
accomplished with the Wang-Landau algorithm. 



We fit the Wang-Landau data (i.e. the second row in Table VII) to the form 



bt{N) = 6,(oo) + ^ + ^ , (5.7) 



and find the fit parameters listed in Table VIII The fits are of reasonable quality, and 



find values for 6t(oo) which lie well below the estimate 6buik — 0.36. It is hard to quote a 
significance for this discrepancy, since we do not have a good estimate of the error in 6buik- 
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FIG. 16: The strong-to-weak transition coupling, bt, plotted versus 1/A^. [Red] squares show our 



results using the Wang-Landau algorithm from Table VII while [blue] circles are from Refs. |1 H I12 1 
[131 [U]. The solid black curve is the fit described in the text. The [magenta] dashed horizontal lines 
give the range in which the bulk transition in SU{12) should take place, according to hysteresis 
scans performed in a lattice theory |3T]. The insert shows a close-up of our new data. 

If we use the error in our results, the significance is between ~ 50 to ~ 200 a. We thus 
think it is very unlikely that bt in the QEK model can be identified with febuik of the SU{oo) 
lattice gauge theory. 

We have compared fits with and without the 1/A^ term, and find that the former is slightly 



preferred, as shown in the Table. It is this fit which is included in Fig. (16). We have also 
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Type of fit &t(oo) 



A 



B xVd.o.f. 



^ = 0,-B/ 0.3142(2) - -7.59(18) 1.45/1 

0.3148(10) -0.037(65) -7.06(97) 1.1/1 

TABLE VIII: The parameters 6^(00), A, and B, obtained from fitting the Wang-Landau data in 



Table VII to the form Eq. (5.7) 



attempted to fit simultaneously to the Wang-Landau results for = 40 and 50 and the 
(more accurate) FS results for = 20 and 30. This fit fails, quite likely because, given the 
very high accuracy obtained with the FS method, we need to include terms of 0{1/N^). 



VI. SUMMARY AND DISCUSSION 



In this paper we have studied the validity of large- reduction for the four dimensional 
quenched Eguchi-Kawai model. This model is a variant of the original Eguchi-Kawai model 
in which the distribution of the eigenvalues of the link matrices is forced to be uniform by 
quenching, while all other degrees of freedom remain dynamical. 

We find that while enforcing a uniform eigenvalue distribution is indeed a necessary 
condition for large- reduction to hold, it is not sufficient. The reason is that quenching 
fixes the eigenvalues only up to permutations that can be performed independently in the 
four directions. These permutations occur dynamically in the model due to fluctuations in 
the unquenched degrees of freedom, and can lead to correlations between the ordering of 
the eigenvalues of the four link matrices. If such correlations occur then we show that the 
arguments of Refs. [U [TUl |T71 [THl [IHl 120] for the validity of the large- A^ quenched reduction 
break down. 

The question then is whether such correlations between link eigenvalues occur. We show 
that they are indeed expected in the weak-coupling regime by minimizing the free energy 
with respect to the ordering of the eigenvalues. This then leads us to perform a detailed 
numerical study of the QEK model with intermediate and strong couplings using Monte- 
Carlo techniques. We find the weak-coupling calculation is indeed a good guide and obtain 
the following evidence for the breakdown of large- A^ reduction in the model: 

• We observe clear evidence for eigenvalue correlations by measuring order parameters 
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that explicitly probe the correlation between the different link matrices along the 
different Euclidean directions. 

• When we compare the plaquette expectation values of the QEK model and of large 
volume lattice gauge theories, we find very large discrepancies that do not go away 
with increasing A^. 

• When we measure the coupling at which a strong-to-weak transition occurs in the QEK 
model, and compare it to the coupling at which the "bulk" transition takes place in 
large- lattice gauge theories in large volumes, we observe a large discrepancy which 
is of order 13%, and very significant statistically. 

We checked that these conclusions are insensitive to the precise form of the quenched eigen- 
value distribution, and to the way we perform the quenched average. We also considered 
values of up to 200 to look for a late onset of behavior, but find none. We conclude 
that the momentum quenched large- iV reduction of SU (N) lattice gauge theories fails in the 
continuum limit. 

We have focused in this paper on the behavior in the weak coupling region, since this 
is where a continuum limit might be taken. Nevertheless, it is also interesting to consider 
the status of reduction in the strong coupling regime. In the strong-coupling expansion 
no eigenvalue correlations appear and so the QEK model is expected to be equivalent to 
the SU{oo) gauge theory for large enough 't Hooft coupling A. It follows that reduction 
is valid until a transition occurs into a phase in which eigenvalue correlations appear. For 
the quenched Eguchi-Kawai this occurs at the strong-to-weak transition. We have checked 
numerically that the eigenvalue correlations do vanish on the strong-coupling side of this 
transition. A similar picture holds for both the EK and TEK models: reduction holds for 
large enough A but is lost below a certain coupling. We stress, however, that this transition 
coupling differs for all three theories (and also differs from the bulk transition coupling for 
SU (oo j). This is just a reflection of the fact that the weak-coupling phases in these theories 
are unrelated. 

These results, together with those of Refs. [271128] concerning the TEK model, mean that, 
currently, only two single-site models are known that can possibly reproduce the properties 
of QCD at large-A^. The first is the "deformed" Eguchi-Kawai (DEK) model, which is 
the single-site example of a class of models proposed very recently in Ref. [33]. In the 
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DEK, the action of the Eguchi-Kawai model is deformed so that breakdown of the Zff 
symmetry is energetically disfavored, and yet at the same time the large- dynamics is not 
modified. Thus in this model the original Eguchi-Kawai proof of reduction is expected to 
remain valid. In preliminary calculations of the DEK model we see that this deformation 
must include terms that decorrelate the gauge fields in different Euclidean directions, and 
this makes a direct connection to our results in the QEK model, where we see that such 
correlation is dynamically preferred. Deformation comes, however, at a cost. Adding all 
possible deformations is likely to be prohibitively expensive, because there are ~ iV^ in four 
dimensions. Whether one can improve this scaling by wise choices of the deformations is a 
subject for future investigation. 

The other single-site candidate is the model obtained by adding 1 < Nf < 4 Majorana 
adjoint quarks, with periodic boundary conditions, to the Eguchi-Kawai action [31j. Here 
the one-loop potential for the link eigenvalues is repulsive if the quark mass in units of 
the lattice spacing, aiat"^, is small enough, and reduction is expected to hold. By taking 
o-iat^ 1 and yet m much larger than the confinement scale A, this single site model 
should reproduce the pure Yang-Mills theory in the IR, and it remains to be seen how such 
a construction compares computationally to that of [33j. For m <^ A this model describes 
the large- A^ limit of QCD with adjoint quarks which, for Nf < 4, is expected to confine and 
to be related to 3-color QCD through the orientifold large-A^ equivalence. For Nf ~ 5 this 
model is expected to be close to conformal, and also of interest. We leave the exploration 
of both these single-site models to future studies. 
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